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We consider a mixture of passive (i.e., Brownian) and active (e.g., bacterial or colloidal swimmers) 
particles, and analyze the stability conditions of either uniformly mixed or phase segregated steady 
states consisting of phases enriched with different types of particles. We show that in sufficiently 
dilute mixtures the system behaves as if it were exposed to two separate heat baths of uneven 
temperatures. It can be described within a second virial approximation neglecting three body and 
higher order collisions. In this approximation, we define non-equilibrium “chemical potentials” 
whose gradients govern diffusion fluxes and a non-equilibrium “osmotic pressure”, which governs 
the mechanical stability of the interface. 


Introduction. Suspensions of actively moving parti¬ 
cles, performing mechanical work at the expense of inter¬ 
nal or external energy consumption, have attracted much 
attention over the last years m- The interest is moti¬ 
vated by biological applications, but these studies shed 
also light on the fundamentals of statistical mechanics. 
These systems share many interesting properties such as 
spontaneous flows uni, but one of the most exciting phe¬ 
nomena is phase segregation It is often driven 

by variants of the so-called quorum sensing, which is the 
feedback mechanism reducing the activity of a given ac¬ 
tive particle in the presence of a high concentration of 
other active particles. Another type of phase segregation 
can occur in a mixture of particles with different levels 
of activity (active and non-active), when one phase is 
enriched in active and the other one in passive particles 
jTTl IT2] . This type of active phase segregation is far less 
understood. 

A new spin on the problem comes from Ref. j!3| . 
These authors study eukaryotic nuclei and the spatial 
segregation between eu- and hetero-chromatin, i.e., be¬ 
tween actively processed and almost silent parts of the 
genome. Viewing the genome as a polymer, they argue 
that genes which are being expressed and, therefore, sub¬ 
ject to RNA polymerization and other active processes, 
should be viewed as active monomers, and silent genes 
are passive monomers. They then hypothesize that the 
observed compartmentalization between the two kinds of 
chromatin is a phase segregation or rather a microphase 
segregation [14] based on activity. Their computational 
model appears reasonably consistent with the data, even 
though the simulation replaces the “activity” by a suffi¬ 
ciently high effective temperature imposed on the active 
monomers by a separate heat bath. One striking obser¬ 
vation is that the effective temperature must be signif¬ 
icantly higher than the real temperature of the passive 
monomers (by about a factor 20). 

In this work we develop a minimal analytical model of 


phase segregation between active and passive particles, 
call them A and B. Our main idea is to look at systems 
of sufficiently low concentration, where we resort to the 
type of reasoning which for equilibriums systems leads 
to a virial expansion. We give a systematic development 
of the second virial approximation in which only pair 
collisions between particles are considered. 

If the interactions are short range, at low concentra¬ 
tion, each particle completely looses its orientational cor¬ 
relations in the time between collisions with other parti¬ 
cles. We can then consider the particles with two distinct 
levels of activity as exposed to two different heat baths, 
with temperatures T4 ^ Tg. This is clearly a system far 
from equilibrium, with energy flowing from the hotter 
to the colder reservoir via the interactions between par¬ 
ticles; more physically, energy is taken from the source 
of activity and dissipated into the surrounding medium 
via our system of particles. There are many examples of 
systems whose description involves two distinct temper¬ 
atures, ranging from plasmas (see basics in, e.g., [15]). 
to spin glasses [IB] and heteropolymers HU. Other ex¬ 
amples are given in [M]. In all these works (with the 
notable exception of Ref. EU), the two temperatures are 
used to describe motions on vastly different time scales. 
In our system, there is no such separation of time scales, 
and our goal is to explore how the access to two different 
heat baths enhances the tendency towards phase segre¬ 
gation. 

Our starting point is the over-damped Langevin equa¬ 
tion 

CiX i = -d i U+(2T i Ci) 1/2 ^(t) , (1) 

for every particle i in the system. Here Xi indicates the 
position of particle i (for brevity, we will make no distinc¬ 
tion between particles in 1 or 3 dimensions). We assume 
that all forces acting on a particle derive from a poten¬ 
tial energy U, while di is the derivative with respect to 
Xi. The friction coefficient for particle i is Q, and £i(t) is 
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a standard zero mean and unit variance Gaussian white 
noise, independent for all particles i. Finally, T) is the 
temperature of the heat bath interacting with particle i, 
it is either T 4 or Tg. 

A prototypical system that may be described by this 
model include mixtures of actively swimming bacteria 
with either oxygen-starved ones (which do not actively 
swim), or just similar sized inactive colloidal particles; 
another example is a mixture of passive colloidal particles 
with the ones capable of light-induced catalytic transfor¬ 
mation of the solvent, like in [4] . In the jargon of the field 
[221 23 , these are either run-and-tumble or active Brow¬ 
nian particles models, sufficiently dilute to be amenable 
to a virial approximation. Another major assumption 
behind our model is that we ignore hydrodynamic inter¬ 
actions (Rouse model), each particle experiencing a local 
friction against an immobile solvent. 

Two particles. We first examine a system with only 
two particles, one A and one B. In this case, Eq.Q 
consists of two coupled Langevin equations. Dotsenko et 
al [2T] studied this problem when the potential energy 
U (a"_ 4 , xb) is a positive-definite quadratic form of its two 
variables. We only assume here that U depends on the 
distance between particles, r = xa~xb, he., U = rr 4 B (r) 
and consider mostly cases where rr 4 e (r) vanishes at large 
r. We derive a Langevin equation for the variable r, by 
combining the two equations 0: 

C,r = F(r) + (2 C r T) lf % , (2) 

where the relevant friction is ( r = (a(b/((a + Cb) and 
the relevant temperature is the mobility-weighted aver¬ 
age T = (CbTa + CaTb)/(Ca + Cb)- The definition of the 
effective temperature T is dictated by the condition that 
the noise £ r (i) is a zero mean, unit variance Gaussian 
white noise. Since C,a,b > 0 are positive, T is always be¬ 
tween Ta and Tg. It follows from the Langevin equation 
([ 2 ]) that the relative distance between particles, r, in a 
steady state is Boltzmann distributed with the average 
temperature T, despite the fact that the system remains, 
of course, out of equilibrium: 

P(r) = exp {-U(r)/T) z” 1 , (3) 

x is here the “partition sum” ensuring normalization m- 

Fokker-Planck equation, currents, and viola¬ 
tion of detailed balance. The Langevin equations 0 
can be recast as a Fokker-Planck equation for the joint 
probability distribution of the coordinates of all particles 
P({a;}) and the corresponding currents <p: 

P = -d z J z , J x = -diU P/0 - TAP/Ci , (4) 

(using the Einstein convention for the summation over 
repeated indices). At steady state, for two particles 
the probability P(r) depends only on the distance r = 
xa — xb- Then 8aP = —<9gP, as well as 8aU = —ObU. 
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FIG. 1: Field of currents. For illustration purposes, the inter¬ 
action potential is chosen in the form u(r) = l/(r 2 + 1), and 
temperatures are Ta = 0.5 and Tg = 1. In the enlarged win¬ 
dow, drift (upward and to the left) and diffusion (downward 
and to the right) are shown separately; unlike equilibrium sys¬ 
tem, these currents are not collinear, their vector sum gives 
rise to the non-potential current field. 


Furthermore, under these conditions, there is no cur¬ 
rent in the r = xa — xb direction. The current vector 
J must be in the perpendicular direction, which means 
Ja~ Jb = 0 , or 

(-d A U/( A + 9 bU/(b) — {T A d A P/ (a - Tb^P/Cb) = 0 . 

This, of course, reproduces the Boltzmann distribution 
with the average temperature (|3|. We now compute the 
current J 4 , taking advantage of the above result to re¬ 
place the diffusion term: 


Ja = Jb 


Ta-Tb 
TaCb + TbCa 


P{r)d A U . 


(5) 


As expected, the current vanishes for an equilibrium sys¬ 
tem (Ta = Tg), where detailed balance is obeyed. How¬ 
ever, if Ta p Tb , detailed balance is violated. For in¬ 
stance, if the system is 1 -dimensional, then its configura¬ 
tion space (xa, xb) is a 2 -dimensional plane, and it is easy 
to visualize loops of current J in such a plane as shown 
on Fig. |T] Physically, as Eq 0 suggests, these loops of 
current mean that the more passive particle moves, on 
average, mainly in the direction of the force acting on 
them, while the more active particle moves in the direc¬ 
tion opposite to the force acting on them. 

Power transfer. As we pointed out above, en¬ 
ergy is transferred everywhere in the system from the 
“hot heat bath” to the “cold” bath, or from the en¬ 
ergy source of active motion to the surrounding passive 
medium. The transferred power from B to A is the aver¬ 
age w = (—x^djJJ)- We note that velocity xa = Ja/P, 
while the average involves an integration with a weight 
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P. Therefore, w = — f J A d A Udx A dx B , leading to: 


WB^A 


T b ~T a 
TaCb + TbCa 


(■ d r u AB (r)Y 


S (r)/T 


-dr . 


. (6) 

As expected, this power transfer from B to A vanishes 
for the equilibrium system if Tg = T 4 , but it is positive 
if Tg > T 4 and negative otherwise. 

To understand the meaning of the result © , consider 
a 3- dimensional system (when the integration over dr 
runs over the volume) and an interaction potential U 
that does not bind the particles together. An example is 
a repulsive U (r) with a little bump of energy scale Uq and 
spatial scale i. Then the integral in Eq. © is estimated 
as (Uq/£) 2 £ 3 / L 3 , with L 3 the box volume (which enters 
in z), and then (assuming for simplicity C-4 = Cb = C) 
the result can be re-arranged as 


—2 


w~[t (T/C) /L 3 } (2g - T A ) Uq/T' 


(7) 


The first factor in the square brackets is the inverse 
Smoluchowski time between collisions of two particles, 
and, therefore, the second factor is an estimate of the 
energy transferred during one collision. 

Many particles. Consider now a system of N, par¬ 
ticles i (■ i,j = A,B). The Fokker-Planck equation 0 is 
generalized for any number of particles. Integrating the 
Fokker-Planck equation over all coordinates except for 
one, we derive a diffusion equation for the single particle 
probability (proportional to the concentration) for every 
particle species: 


dtf{r) = Na 
dt Ca 


du AA 

dr 


pf A 


(r, r') dr' 


+ 



du AB 

dr 


pi B 


(r, r') dr' 


Ta 

Ca 


Vfrftr) ■ 


( 8 ) 


A similar equation is obtained for p B (r) . These equations 
for single particle probabilities include pair probabilities 
P 2 ■ By integrating the multiparticle Fokker-Planck equa¬ 
tion over all coordinates except for two, we derive equa¬ 
tions for P 2 , which include the 3-body correlations P 3 , 
and then a hierarchy of equations |28| . However, if the 
density is small enough, we can neglect the 3-body corre¬ 
lations and ignore all terms involving p 3 , thus obtaining 
closed equation for p 2 . Consistent with the two particle 
system, we obtain in this approximation: 

P 2 ( r > r ') = P\{*)P J i(r') exp [-u lJ (r - r') /fa] (9) 


where i,j = A,B. The effective temperatures entering 
these expressions is different for the three types of inter¬ 
actions: T aa = T a , T B b = Tg and T AB = T. 

It is important to note that these distributions form 
only as a result of an averaging over many collisions hap¬ 
pening in the system under steady state conditions (sim¬ 
ilar in this respect to an equilibrium system). 


Inserting the ansatz (34) into equations 


and in¬ 


troducing the concentrations c*(r) = Nip\( r), we obtain 
closed equations for the concentrations: 


dc l (r) _ / jdpA 

dt Ci dr \ dr ) 


( 10 ) 


These equations look like regular diffusion equations, but 
they are governed by non-equilibrium analogs of chemical 
potentials [29l : 

PA = T a In c - 4 + T a B a c a + TB ab c b (11) 


and a similar equation for /xg. The virial coefficients 
are defined each with its own temperature, as B^ = 


/ 


1 _ e -“ u ( r )/ r .j 


d 3 r. 


Non-equilibrium chemical potentials, as quantities 
whose gradient determines the flux, were discussed in 
Ref. [5:. It was shown, that, unlike its equilibrium coun¬ 
terpart, a non-equilibrium chemical potential, in gen¬ 
eral, cannot be obtained as a derivative of a free en¬ 
ergy. This was shown in particular for the gradient terms. 
In our case, the situation is different, because the non¬ 
equilibrium chemical potentials p A and /ig appear to be 
the partial derivatives pi = ^ of a function, which looks 
like a two-temperature free energy (per unit volume): 


/ = T A c A In (c A /e) + T B c B In (c B /e) + 

+ (1/2 )T A B A c\ + (l/2)TgHgc| + TB ab c a c b . 


Instability of the uniform state and “spinodal.” 

Suppose that Cq and c B are the averaged spatially uni¬ 
form concentrations of both components. By introducing 
small space dependent perturbations c*(r) = Cq + Sc l ( r), 
we perform a linear stability analysis in the standard way. 
This shows that an instability occurs macroscopically un¬ 
der the condition 


&A T a T b B 2 ab 


1 + < f> A 1 +1 


> 


T B a B b 


(13) 


where we have defined the volume fractions c l Bi = fa. 
In general, the virial coefficients depend on temperature 
in a complex way. A simple limit corresponds to purely 
excluded volume interaction potentials such that the B's 
do not depend on temperature. We study this case in the 
following. 

In the plane (f> A and 4> g, the non-equilibrium equivalent 
of the spinodal line (13) is a hyperbola (see figure [ 2 ]) . The 
contrast between temperatures favors instability, it works 
in the same direction as contrast between interactions. 
But this instability, to have a physical meaning, must 
occur at <p A < 1 and </>g < 1 ; moreover, <p A + <p B < 1 . 
For instance, consider the most symmetric case of iden¬ 
tical particles in all respects except driven by different 
temperatures: Gl = Cb, Ba = B b = B ab . In this case 
the spinodal line is in the physical range (f> A + <t >g < 1 
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<!>A ~ <f>B 



FIG. 2: Phase diagram, <p A and <f>B are volume fractions of 
A- and Z3-particles, while 1 — </>,4 — cj>B (which is the distance 
to lower side of the triangle) is the fraction of solvent. Green 
line is the “spinodal” (below this line, uniformly mixed state is 
unstable). Blue line corresponds to constant osmotic pressure, 
it shows the possibility of two states coexisting. 


if the ratio of the two temperatures is outside the range 
17 - 12 V 2 « 0.029 < < 17 + 12^2 « 34 [3D]. Thus, 

a numerically large temperature contrast is required to 
achieve an instability by temperature difference alone. 
This is somewhat consistent with numerical observation 
of Ganai et al HE as they used T B /T A = 20. 

Pressure and “binodal.” To address not only 
the loss of stability of the uniform mixed state, but 
also the steady state phase segregation, in addition to 
non-equilibrium chemical potentials we also need a non- 
equilibrium equivalent of the osmotic pressure. Given 
that non-equilibrium chemical potentials are the deriva¬ 


tives of a “quasi-free-energy” ( 12 ), we can expect the 


osmotic pressure to be given by Gibbs-Duhem formula 
P = c a p-a + c B /rg - /, yielding 


P = T a c a + T b c b + 


(1 /2)T A B A c z A + (1/2 )T b B b 4 + TB AB c A d 


,AJ3 


(14) 


We have directly derived this result in two independent 
ways, first by computing the force exerted on the wall by 
replacing the wall with a potential ramp, not necessarily 
the same for both particle species, uj^r) and u B r)), and 
second by relating the pressure to the pair correlation 
function in the bulk [241 125] . Since we do not consider 
non- spherical particles which experience a torque upon 
interactions with wall, we do not have the complications 
studied in the recent works [ 22 J [23| t2B] , and, indeed, the 
two derivations |31j yield identical results for the osmotic 
pressure (14). This result means that the densities of 


two coexisting phases at steady state are found by the 
Maxwell common tangent construction based upon the 


quasi-free energy function (12). So far, we make this 


statement in the second virial approximation only. It 
remains to be seen whether this is still true or not when 
higher order collisions are taken into account, which is 
necessary at higher densities. 


The calculation of the dissipation @, can be general¬ 
ized to the dissipation per unit volume of the solution. 
The result is 


w = c A c B 


f du AB y . 


Tb~T a 
T A (b + TbCa J \ dr J 


5 (r) 


d 3 r . 


(15) 

It shows that in a phase separated system the dissipation 
mostly happens around the phase boundary |32j. 

Conclusion and discussion. To conclude, we first 
have to estimate the temperatures T A and T B in terms 
of real parameters of active particles. Since we as¬ 
sumed that the re-orientation time of one particle is much 
smaller than the time between collisions: T r <C r c , the 
particle trajectory between collisions is that of a ran¬ 
dom walk, characterized by an effective diffusion con¬ 
stant D e ff = T e fi/( ~ UqTv /6 + D, where vq is the 
swimming speed, D = T/C, is the passive diffusion con¬ 
stant, T = T b is the real ambient temperature, and 
T e g = T a corresponds in our theory to the tempera¬ 
ture of the hotter heat bath. The temperature difference, 


which controls dissipation rate (151, is therefore given by 
T a — Tg ~ r!gT r ^/ 6 , it is directly related to the level of 
activity measured by the swimming speed vq. 

Returning to the applicability condition r r <C t c , for 
particles of size b with typical distance between particles 
d, the collision time is estimated by the Smoluchowski 
formula r c ~ D e sb/d 3 ~ D 0 ff</)/ 6 2 , where <j) ~ b 3 /d 3 is 
the volume fraction of particles. The condition of appli¬ 
cability is then conveniently formulated in terms of the 
Peclet number vq T r /b = Pe <C 1 />/<j>. Thus, our theory 
should work if the system is dilute enough and/or the 
active drive is not too strong. 

Another significant limitation of our approach is the 
fact that we neglect hydrodynamic interactions. This 
might be particularly important in the case of actively 
swimming colloids, as they usually drive themselves by 
creating and maintaining a train of diffusing chemicals, 
and accordingly their interaction upon approach and col¬ 
lision is hardly describable in terms of a conservative force 
potential, as we did here. 

Despite all limitations, we believe that the theory de¬ 
veloped here is useful because it is physically transparent 
and may be instructive as a source of physical intuition 
for these highly unusual driven systems. 
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Activity induced phase separation: Supplementary material 


In this Supplementary Material, we provide the following: 

1. For two particles, we analyze their joint diffusion and show how it is enhanced by the uneven driving tempera¬ 
tures. 

2. For many particles, we provide a detailed derivation of the hierarchy of the equations for the correlation functions 
and show how our second virial approximation comes out. 

3. We provide detailed analysis of the instability conditions when not only temperatures are uneven, but so also 
frictions and interactions. 

4. We provide two methods to derive the osmotic pressure. 

5. We derive a more general expression for the transfer of power between the heat reservoirs that includes two-and 
three-body correlations. 


Two particles: “Center of friction” diffusion . 

In the main text, we analyzed the relative motion of 
two particles A and B by looking at the variable r = 
x A ~ x B . It is interesting to find also how the presence 
of two distinct temperatures affects their joint diffusion 
in space. To do so, it is convenient to define their joint 
coordinate R in such a way that the noises in r and R are 
statistically independent. This is achieved by choosing 


here the power the spectrum of the force F(r) at zero 
frequency is 

/ OO 

(F(r(t))F(r(t + r)))dr . (19) 

-OO 

in order to find the power spectrum of the force, we 
Fourier transform the Langevin equation for r (Eq.(2) of 
the main text): 


R = 


CaTb 


-xa 


CbTa 


CaTb + C bTa Ca T B + CbTa 

Then, the Langevin equation for R reads 


xb 


(16) 


((A+(B)R= J ^ 1 ^F(r)+ 


+ \/2 (Ca + Cb) ^ ^ 8 £r(0 


(17) 


where £_r(£) is a zero mean unit variance Gaussian 
white noise (independent of £ r , as stated), and T = 
( CaTb + CbTa) / (Ca + Cb) was defined in the main text. 
In an equilibrium system, at X ' A = Tg, the inter-particle 
force F(r) does not couple to the joint motion. Not so out 
of equilibrium: since on average the force F(r) vanishes 
(by symmetry), ( F(r )) = 0, the F(r) term provides an 
additional noise driving the diffusion of the variable R. 
Of course, it is not a white noise, so that the dynamics of 
R is not a simple diffusion on time scales shorter or com¬ 
parable to the correlation time of r. But on longer time 
scales R undergoes simple diffusion, with a diffusion co¬ 
efficient which can be directly read out of the Langevin 
equation (Jl7|) , because F(r(t )) is statistically indepen¬ 
dent from £#(£): 


= T a T b 1/ T a -T b V (f 2\ 

CaTb + CbTa 2 \CaTb+CbT a ) K '“=° ’ 

(18) 


F u = iu( r r u - (2C r T) 1/2 , (20) 

Then multiplying it by the complex conjugate and as¬ 
suming [ii)Fj w _ 0 = 0, we obtain (F 2 )^ =0 = 2 ( r T. This 
yields 


D T A T B 1 (T A -T B f ( aCb 

R CaTb + CbT a 2 CaTb + CbTa (( a + (b) 2 

( 21 ) 

We see that the difference in temperatures, indepen¬ 
dently of the sign, enhances the joint diffusion. 

It is also instructive to write the Fokker-Planck equa¬ 
tion in terms of the variables r and R: 


dtP = 


Ca + Cb 
. CaCb . 


d r ( F(r)P ) - 


T b -T a 


CaTb + CbTa 
CaT b + CbTa 
CaCb 


F{r)d R P+ 

d 2 r P 


t a t b 

CaT b + CbTa 


9rP ■ 


( 22 ) 


Here, the right hand side has explicitly the form of a 
divergence, and can be written as 


d t P = -d r J r - d R J R , 


(23) 
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where the components of the flux are 


Define single particle probability densities, two particle 
probability densities, etc as: 


Jr = 


Ca + C B 

CaCb 
Tb~T a 
UTb + C bTa 


F(r)P - 

F(r)P - 


CaTb + CbTa 
CaCb 

T a T b 

XaTb + CbT a \ 


d r P , (24a) 

Or P ■ 
(24b) 


Hierarchy of equations for the correlation functions 


PiX) = J 5(rf~r)P ({r}) d{r} (27a) 

Pi W = J S (rf - r) P ({r}) d{r} (27b) 

P2 A ( r, r 0 = j 6 ( r f ~ r ) S ( r f ~ r 0 P (M) d i r I 

(27c) 

P BB ( r, r ') = J S ( r ? - r ) 6 ( r f - r ') p (M) 4 r } 

(27d) 

pf B ( r, r') = J 5 (rf - r) 6 (rf - r') P ({r}) d{ r} 

(27e) 


Notations, definitions and symmetries 


As stated in the main text, we operate with a multi¬ 
dimensional Fokker-Planck equation for the probability 
density P as a function of the positions of all the parti¬ 
cles in the system. More specifically, consider a system of 
N a particles A and Nb particles B ; their coordinates are 
rf, rf, , rf A and rf, rf ,..., rf, B . The potential en¬ 
ergy includes single particle potentials and pairwise ad¬ 
ditive interactions: 


U(rf, 


A „A 


r.A , r B r B 


■2 i • • • > 1 Af A i 1 1 ’ 1 2 j • • • ) 1 N b I 
Na N b 

= Y uA ( r i 4 ) + Y u i ( r : 


i 


Na 


2 Y uAA ( r ^ - r ^) + 

Na N b 

YY uAB ( r f~ rB ) + 

i j 
Nb 


1 


Y u BB (rf - r 


(25) 




The probability density P is also a function of all the co¬ 
ordinates: P = P(rf,r A , 


is normalized: 


r A . r B r B 

[ Na. Nt3 ) 


There are 4 types of 3-particle densities: p AAA ( r,r',r"), 
pf AB (r, r', r"), pf BB ( r, r', r"), pf BB ( r, r', r"). All these 
densities are independent of i and j etc, i.e., the prob¬ 
ability density is the same for every particle of a given 
species. 

There are several normalization and symmetry prop¬ 
erties (A, y = A, B): 


/ Pi(r)dr = 1 , 

(28a) 

* 

J-T 

II 

* 

'•-T 

(28b) 

, 10 * 

II 

* 

(28c) 

j pf y {r,r')dxdx' = 1 , 

(28d) 

1 P$ y ( r,r')dr' = pf ( r) , 

(28e) 

[ P* y (r,r')dr = pf (r') 

(28f) 


Fokker-Planck equations for the densities 

Integrating out all variables except for one, or except 
for two, etc, we obtain the following dynamic equations 
for the densities: 


dpf(r) 

dt 


gv^(r) + i-S r 


duf( r) A 


dr 


pf( r) 


N a - 1 
(a 
Nb a 


dr 


du AA (r, r') aa 


dr 

du AB ( r,r') AB 


pf A (r, r') dr' 


dr 


P A ( r , r 0 dr ' 


| + 

(29) 



r A 

C 2 5 ■ 


r (v A ;r?,rf,...,rf, e )d{r} = 1 . (26) 


A similar equation is obtained for p f (r), which is not 
given for brevity. For large numbers of particles, we can 



replace TV 4 — 1 ~ TV 

^^ = g(v? + vf.)^( r,r') + 


■£* 

■g* 


—^P 2 (r.r') 


au^\ A, (r y } 

■/ a,yy, r " ) ^ B ( rr '.r„ )(ir „] + (30) 


chi 7 (r') 


dr 


—p^(r,r ) 


du AA (r,r') AA . 

dr' P2 (r ’ r) 


AU 

Ca 

TVg 

C -4 


d r ' 

d r ' 


7 a “77' r "W "V,r>")ar 

jojq±n p AA B{ryy , )dI , 


A similar equation (not given here) is obtained for p BB . 
But the equation for the mixed probability deserves to 
be written down: 

v = f gi v ; + ?;d ^(,, ,->+ 


dt 

is 

is 

Na 
C A 
TVg 

a 
> 


Ca 

^(r.e 

————————— v AB (r r') 

dr p 2 (r ’ r) 


Cb 


+ 


' gn- 46 (t,t") i _ab Bi 


dr 

duf (r') AB r r / ) 
dr' P2 (r ’ rj 
(r,r') . 

dr' P2 (r,r) 

-joA^y) pSBA[tyyt)dt 


V SB (r,r',r")dr" 


(31) 


AU 

Cb 

TVg 

Cb 


dr- 

dr-' 


The equations for £>3 involve P 4 , and so on, ad infinitum. 

If the density is small enough, we can neglect all triple 
collisions, i.e., directly discard all terms involving p 3 . In¬ 
deed, given the normalization, any term containing P 2 in 
equations (301 or (301 is of order 1/V 2 , while every term 
containing p 3 is of order N/V 3 . This, of course, simpli¬ 
fies the equations quite dramatically, and reduces them 


essentially to what we obtained for the case of two par¬ 
ticles. Omitting the single particle potential terms, we 
obtain 


T?^7 , A> = |(v; + v;,)p#y r ,r') + 


^ 3 - 

i~S 


du AA { r,r') AA 


dr 


P AA ( r , r') 


(32) 


du AA (r, r') AA ' 
-- -p AA (r,r ) 


dr' 


(and a similar equation for p BB ), 
d P f B (r,r') _ /T4 2 T« 2 


dt 


C A V2 r + ^V 2 r ,} pAB (r,r')+ 


1 

(a 

1 

Cb 


du AB (r, r') AB 


dr 


Pf(r,r') 


(33) 


3 “"’ (rir 'W»(ry> 


dr' 


These equations are simple enough to guess the solution 
based on our knowledge of the two-particle case: 


P AA ( r i r ') =p7(r)p7( r, )exp 


pf B (r,r') = Pi (r)pf (r') exp 
P AB { r, r/ ) = P A { r )Pi ( r 0 exp 


..AA 


(r — r') 


,BB 


Ta 

(r — r') 


u AB (r — r') 


T 


(34a) 

(34b) 

(34c) 


Of course, the central feature of this result is the appear¬ 
ance of the average temperature, as defined in the main 
text 


T = 


CaTb + Cb Ta 
C A + Cb 


(35) 


In order to obtain the result (34), we look for a solution 
of the form p xy {r,r') = q xy (r, r') exp [— /3u xy (r, r')], 
plug it into equations (32) or (33), and discover, that 


in the remaining equation for q , the variables separate, 
meaning that q factorizes into a factor that depends only 
on r and a factor that depends only on r'. 


Diffusion equations and non-equilibrium chemical 
potentials 


Plugging the ansatz (34) into equations (291, we obtain 


closed results for the densities. This involves the integral 
f du AA (r,r') , . .. "- 4 - 4 (---') 


A i„'\ 


dr 


-KV) 


A dr' , (36) 


which can be integrated by parts. Finally, we obtain 
equations that look like diffusion equations for a regular 
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system in contact with a thermostat, 

0g*(r) l d ( A dp A \ 

dt ( A dr \ Or ) 
dc B {r) _ 1 d f B dp, B \ 
dt Cb dr \ dr ) 


(37a) 

(37b) 


but these equations contain non-equilibrium chemical po¬ 
tentials ■, as stated in the main text (Eq.(ll)). 


Linear stability analysis 

Suppose that c A and c B are the averaged spatially uni¬ 
form concentrations of both components. By introducing 
small space dependent perturbations c- 4 (r) = c A +Sc A (r) 
and c e (r) = c B + Sc B ( r), we perform a linear stability 
analysis in the standard way: 



FIG. 3: Spinodal line ( |39| ) (in green), shown for one particular 
choice of parameters (Ta/Tb, (a/(b, BaB b /B ab )- Physical 
meaning has only region (j> A + (j>B < 1 > the instability region 
within this region is shaded. 


{^r = C 7 y2 [( Ta + t a4Ba) 5c a + ( t c a b ab ) <5c B ] 
1 nf = h v2 1(Tc b B ab ) 5 c a + (T B + T b c b B b ) <5c e ] 

(38) 

This shows that an instability occurs (at q = 0, i.e., 
macroscopically) under the condition that the determi¬ 
nant of this matrix vanishes, i.e., the system is unstable 
if 

c a c b T 2 B 2 ab > T A T B ( 1 + c a B a ){ 1 + c b B b ) . (39) 


At the instability, the unstable combination (eigenvec¬ 
tor whose eigenvalue flips sign) is 




5c B (r)C, B 


T A B A c A (1 + B a c£) JT b B b c b (1 + B b c b ) 


(40) 


In the plane (c A , c B ), the spinodal line (391 is a hy¬ 
perbola (Fig. [3]). A better way to represent it is to use 
a triangular phase diagram as given in Fig.2 of the main 
text. 



FIG. 4: Instability exists above this surface, see condition 
(441. 


• Contrast of frictions, 


Non-equilibrium “spinodal” line for athermal 
particles 


_ Ca - C B 
(a + C B 


-1 < K < 1 . 


(43) 


The dimensionless parameters of the system are as fol¬ 
lows: 


The condition that the spinodal exists within the phys¬ 
ical range <p A + (j> B < 1 reads 


• Contrast of excluded volumes, 


P 


B AB 

\Jb a b b 


• Contrast of temperatures, 


t a -t b 
t a + t b 


— 1 < T < 1 . 


(41) 


(42) 


4 k 2 t 2 — 4 KT + 1 

1 — T 2 


/3 2 >9 . 


(44) 


This condition is presented graphically in two different 
ways, in Fig. [4] in the form of a 3D surface and in Fig. 
[5] as an array of 2D plots. Beautifully, the contrast of 
frictions n becomes irrelevant for the equilibrium system, 
when r = 0 . 
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FIG. 5: For every k, instability exists above the line, accord¬ 
ing to Eqns (441. 


Non-equilibrium osmotic pressure and “binodal” 

To address the steady state phase segregation, in ad¬ 
dition to the non-equilibrium chemical potentials we also 
need to define a non-quilibrium osmotic pressure. We 
derive it in two different ways. 


Derivation 1 


To find the osmotic pressure, imagine that the system 
“feels” single particle potentials u A (r) and uf(r) such 
that they are both like a box, except that one wall of this 
box has a (not necessarily very) sharp potential “ramp” 
in the direction, perpendicular to the wall: u A,B (r) = 
f A , B x, as shown in Figjb] In this case, the pressure is 
found according to 


POO POO 

P = fA / c A {x)dx + f B c B (x)dz 
Jo Jo 


(45) 


because every particle A present in the ramp area exerts 
on the wall the force f A , and similarly for B. We em¬ 
phasize, that this is actually an osmotic pressure, in the 
sense that the ramp potentials u A,B (r) = f A , B x act only 
on the A and B particles while the solvent penetrates 
everywhere completely freely. This means that our ramp 
potentials represent an osmotic piston. 

A similar expression for the osmotic pressure was also 
used in Ref. S5|, where it is derived from the expression 
of the Helmhotz partition sum, i.e., from equilibrium sta¬ 
tistical mechanics. We feel necessary to emphasize that 


Eq. (451 is derived on purely mechanical grounds, and it 


has nothing to do with thermodynamic equilibrium. As 
such, it is perfectly applicable to our present problem. 

To find the steady state concentration profile in the 
presence of ramp potentials, we slightly generalize the 


diffusion equations (371 by including the external poten- 


At steady state, the concentration profile must be such 
that /z + u\ = const for both the A and B components. 
This can be written as 


c A (r) = C A e 


1 + c a (v)B a + c b {v)B ab 


T 

T~a 


c B (r) = C B e t b 


1 + c B (r)B B + c a {y)B ab 


T 

Tb 


(47a) 

(47b) 


Here C A and C B are normalization factors. To make 
things simple, we assume that the “ramps” are not too 
shallow, such that the normalization integral is domi¬ 
nated by the bulk volume V where both ramp poten¬ 
tials vanish. Given that the virial terms in the chemical 
potentials are the corrections to the ideal gas, we solve 
iteratively and get 


l (r) = 


Na _El. 
— e 


N B „ T 
+ B ab — 
V t a 


AU 

V 


B A l-e~ 


1 — e 


(48) 


and similarly for c B ( r). Note that the result does not 
depend on the ramp forces f A and /g, which do not have 
to be identical. 


Derivation 2 

Our starting point of the second derivation is the ki¬ 
netic expression of the pressure 

P = Pideai - ^7 ( r b' ' d nj u ( r a) / • ( 49 ) 

Sometimes it is called Irving-Kirkwood formula (21, ]25]. 
As in the first derivation, the important point is that this 
equation follows from pure mechanics, does not make any 
assumption related to equilibrium statistical mechanics. 
In terms of pair distributions the Irving-Kirkwood 
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formula reads 


particles: 


P =Pideai - ^ j r 12 ■ d ri2 u AA (r 12 )p^ A (r 1 ,r 2 )~ 

_ N a N b j Fl2 . dri2U AB {ri2)p AB {ri r2) _ 


~~yJ r i2 ■ dr 12 u BB (r 12 )p^ B (r 1 ,r 2 ) . 


(50) 


Using the ansatz 34 for p 2l and integrating by parts (and 
remembering that V • r = 3), we obtain the same result 
as before for the osmotic pressure. 


Power transfer 

In the main text, we outlined the derivation of the 
power transfer in the cases of either two particles, or 
many particles with only pairwise collisions. Here we 
establish a more general result which suggests that the 
power transfer is expressed in terms of only pair and 
triple correlation functions (but not higher order ones). 
Consider the work performed by all forces per unit time 
on all A particles, which is also the power received by A 


Na 


r du 

T a dP ( dU P' 

1 dr A 

[Udr A+ dr A U\ 


d{ r} , (51) 


with U = [/({r^},{rf}) the total potential energy of 
the system. In the integral ([51] ), the first factor is the 
force which acts on particle Ai due to all other particles, 
while the second factor (in square brackets) is the cur¬ 
rent, i.e., the velocity of the particle Ai multiplied by the 


probability density P. Thus, the integral (511 is the av¬ 


erage power transfer to one particle A r . By symmetry, it 
is independent of i, so that the summation over i reduces 
to a factor Na- As long as the interaction potentials are 
pairwise additive (251, the force is also a sum: 


dU 

dr A 


Na q u AA (, C A _ r Aj N a du AB ^A _ r B) 


I A* 


dr A 


n b 

£ 

k 


dr A 


(52) 

The first term is the force acting on particle A t due to 
other A particles, by symmetry this term vanishes on av¬ 
erage in the sum over A particles. And the second term, 
which is the due to B particles on A particles, yields: 


W= N^ j aii- 4 ^ (r - r') + ^ /(AdiUAT) ^ {r y )Ml ' + 

(a J dr dr (a J V dr J 

Nj^ r du AB (r — r') du AA (r — r") 


C A 


dr 


dr 


-p ABA (r,r',r")d 3 rd 3 r , d 3 r"+ 


(a 


dr 


dr 


(53) 


NaNb f du AB (r — r') du AB (r — r") 

- X — 


n ABB 
P 3 


(r, r', r") d 3 rd 3 r'd 3 r" 


Neglecting the three body collisions (terms with p 3 ), and Here, we emphasize once again that, as long as in¬ 
using the known expression (34) for p 2 , we return to the teraction potentials are pairwise additive, as in Eq.(25), 
result given in the main text. 


Eq.( 53) is exact. 































